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Abstract 

The double ratio {fB s /fB d )/(fD 3 /fD d ) is calculated by QCD sum rules in heavy quark 
effective theory (HQET), both numerically and analytically. Our expression for the double 
ratio shows explicitly the dependence on the light quark masses, the heavy quark symmetry 
breaking and the vacuum condensates. The numerical result favors the double ratio to be 
a little greater than 1. 
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Precision determination of the double ratio (fB 3 /fB d )/{fD a /fD d ) is of important interest in 
understanding the B — B mixing phenomena. Both lattice and QCD sum rule calculations [0, 
have given that JbJ fB d and Jd s I fo d ar e close to 1. Grinstein0, based on the phenomenolog- 
ical Lagrangian incorporating heavy quark and chiral symmetries, has found the double ratio 
ifB s /fB d )/(fD d /fD d ) = 0.967. Oakes|4j] has also done similar work by making use of the Wigner- 
Eackart theorem concluding that /bJ fB d = 0.989, /dJ fo d = 0.985 and (fB 3 /fB d )/(fD d /fD d ) = 
1.004. Also, It is noted that there is a difference in the sign of correction term away from 1 
between Grinstein's and Oakes's expressions for the double ratio. Combination of QCD sum rule 
method with HQET, as a powerful approach to nonperturbative dynamics in the weak decays 
of heavy hadrons, provides us with a possibility to discuss the double ratio from QCD theory. 
In the previous paper ||, we have estimated the decay constants of heavy mesons and their 
SU(3) breaking effects by using QCD sum rule in HQET, obtaining fsj ' fs d = 1-17 ± 0.03 and 
fD s l ' fo d = 113 ± 0.03. Therefore, the double ratio can be extracted from the above results as 
(fB a / fB d )/{fD a / fD d ) = 1.035 ± 0.06. Obviously, the indirect calculation brings about a large 
error. To get a more precise result, in the present work we try to discuss it directly within the 
QCD sum rule framework of HQET. It should be emphasized that in our approach the SU(3) 
symmetry breaking effect in the double ratio, as has been shown in Ref. 0, depends on not only 
the current quark masses, but also the condensates < 0|gg|0 > and < 0\qaGq\0 >. This is 
different from the investigations of Grinstein and Oakes, where the SU(3) symmetry breaking is 
determined only by the current quark masses. 

let's begin with a brief review on Grinstein's and Oakes's works. Decay constants of heavy 
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mesons M are denned as 

< 0\q lfll5 Q\M >= if M P^. (1) 

According to Grinstein||, the ratio of decay constant, for instance, fBs/fBd , can be expanded in 
a quantity m s /A(A is QCD scale) if only SU{2>) breaking is taken into account, fsj fB d = 1 + ^- 
Heavy quark symmetry breaking modifies this expression as fsj fB d = 1 + ^( a o + a i^" + ' ' ')• 
The double ratio incorporating SU(3) and heavy quark symmetry breaking effects may read 
approximately 

UbJ '/bJ/C/W '/dJ = 1 + ai(m s /m 6 - m s /m c ). (2) 

Setting a\ = 1 for the coefficient of the correction term and m s = 150MeV, m c = 1.5GeV and 
mj = 4.5GeV for the current quark masses, the double ratio is approximately 0.93, which is 
smaller than 1. To better estimate the double ratio, Grinstein||, making use of the phenomeno- 
logical Lagrangian for heavy mesons and light pseudoscalars incorporating heavy quark and chiral 
symmetries, found it to be 0.967. In Ref.f|], Oakes, utilizing the relation d^A^ = i{mQ+m q )q^Q , 
obtained exactly 

/a//Bi .f^ ! f^«i (3) 

and 

\M D J \m c + m d J <0\dj 5 c\D d > 

Thus, the ratios of decay constants can be precisely estimated with the help of the Wigner- 
Eackart theorem which allows the approximation relation < 0|s75<5|Mq s >^< 0\dj 5 Q\MQd > to 
hold, to very high accuracy, 



M B J \m b + m d J <0\d'j 5 b\B d > 



SdJSd{ = (^Y(T^) <<^P-> (4) 
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h,/h^l A M(^^) (5) 



M B J \m b + m d 



r ■ r M Dd \ 2 (m c + m s . 

f D Jf Dd = hnr* L. , ^ • ( 6 ) 



and 

fr, /fr. =( 

MdJ \m c + m d , 

From Eq.(5) and Eq.(6), the double ratio (fB s /fB d )/{fD a /fD d ) = 1-004, and an approximate 
expression can be found 

(f If )/(f If 1-1 n( M Bs~ M B d M Ds -M Dd \ m s -m d m s -m d 

To the extent that M Bs — M Bd ~ M Ds — M Dd ~ m s — m rf ~ m s , M B ~ m b and M D ~ m c , Eq.(7) 
is simplified as following 

(WW/(WM)-l + a(^-^), (8) 

which shows that the double ratio is, indeed, unity up to a correction of order a{m s lm b — m s jm^)^ 
albeit with coefficient a of -1 rather than 1. 

Now, we turn to calculation of an analytic form of the double ratio by means of QCD sum 
rules in HQET. The heavy quark effective lagrangian to next-leading order in 1/mq is given 
by§ 

L = h^-Dh v + — + (9) 

where L K = h v {iD) 2 h v is the kinetic energy operator of heavy quark , L s = ^hvgsa^G^hv 
is the chromomagnetic interaction operator and C mas (/i) = [£^i^fil] 3 /( 11 - 2n //3) \^ _|_ is the 

Wilson coefficient. To calculate decay constants of heavy mesons by HQET approach, matching 
current operators in fully theory onto ones in HQET is necessary. Let's calculate the axial current 
matrix element between the vacuum and a heavy meson state. A match|| |7| of the axial current 
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operator in fully theory onto the local operators of the effective theory gives the expression for 
the decay constant at order 1/rriQ in the heavy quark expansion in the case of light quark mass 

SuJnS = F(m Q )[l + d M ^]{l + ^> + ^[G,(m s ) - |] } , (10) 

with d,M — 3 for pseudoscalar mesons and —1 for vector mesons. F(mQ),Gj(:( m Q) an d G^(m,Q) 
are some renormalization group invariant quantities^]: 

F(7n Q ) = [C 1 <ji) + ^]F(ri, (11) 

A A 777 

G K (m Q ) = GM - -b{fi) + -±a{fJL) (12) 

O 

and 

A A 777 

G s (m Q ) = G mag (fi)GM - -[£(//) - 1] + ^(/i). (13) 

In the above equations, the short— distance coefficients Cx(ji), C^/x), ^4 (/•*), &(a0 and -B(^) 
have been obtained in Ref.||. F(fj,), Gk(^), Ge(^) and A are some long— distance parameters [|, 
0| : F(p) is used to parametrize the leading order matrix element in I/tuq expansion, Gk(/J>) and 
Gs(/i) are two additional universal parameters, which are respectively introduced to parametrize 
the 1 / rriQ corrections to the hadronic wave function by insertions of the subleading operators L K 
and Ls in the effective lagrangian into the matrix element of the leading order current, and A 
expresses mass difference of the heavy meson and the heavy quark inside it. It is interesting to 
note there is an explicit dependence of these low energy parameters on not only the light quark 
masses but also the condensate parameters in our approach, as shown in Ref.||. Considering 
these effects, we can obtain easily the expressions for the ratios of decay constants 

fBjfB * - F d (, )V ff^ 1 + 9MM (14) 



and 

foJfD * ~ F MV m; 1 + £^ (15) 

TYL C 

with C Q M = G q K (n) + 6C^ as (//)G|(/i) + ^f^m q - b -^f^A q . Following the standard method 
for QCD sum rule calculation]^], one can get the sum rules for those long— distance coefficients 0: 

A, = —^InK, (16) 

FMe~ 2A « /T = K, (17) 



1 d /TJ 



= 55^) (18) 



and 

AdT\ K 

Here, we have not given the final expressions [[J for K, Jk and Js, which are dependent on the 
Borel parameter T, threshold uj c of the continuum states and light flavors. In the light of Eq.(14) 
and Eq.(15) and by a simple manipulation, we have the double ratio 

UBjfB d )(f D Jf Dd ) = 1- G(/i,m b )^ + G(/i,m c )« 

+ AG^m^-AG^m^ (20) 

to next-leading order in 1/rriQ, where £?(//, wig) is defined as 

G{fi,m Q ) = -[B(n,m Q ) - A(fi,m Q )] + ^[b{fi,m Q ) - a(/i,m Q )] (21) 

and AG(fi, tjiq) is SU(3) symmetry breaking parameter 

AG{/i, m Q ) = AG K (n) + 6C ma9 (/i, m )AG E (/i) (22) 

with 

AGM = G s M-G d M (23) 



and 

AGM=G s M~G d M. (24) 

In deriving Eq.(20), the approximate relations M Bs — M Bd ~ M Ds — M Di ~ m s — m d and 
A s — A d — m s — m d have been used. 

The Wilson coefficients G(/i,mq) are easily evaluated via taking \x = lGeV, m b = 4.5 ~ 
5.0 GeV and m c = 1.2 ~ 1.5 GeV. Now we are in a position to compute two low energy 
parameters AGk(^) and AGs(/i), by resorting to the QCD sum rules in HQET. Before numerical 
manipulation, we would like to deal simply with these two quantities. If we assume both G 9 K (fi) 
and to be continuous functions in the light quark mass m q , AGk{h) and AG s (/i) can be 

expanded in terms of {m s — m d ) respectively 

AG K (fi) = f K (m s -m d ) + --- (25) 

and 

AGM = fe(m s -m d ) + ---, (26) 

where f K =^\m q =m d , h = S^\m q =m d and " • • • " denotes the higher power terms of (m s - m d ). 
If we neglect contributions from the higher power terms, Eq.(20) becomes 

(hJhMh.lh.) = 1 + 4^) - 4^) W 

with ai = f K + GC mag (fi, m b )f-E - G(/i, m b ) and a 2 = f K + 6C mag (fi, m c )/ s - G(/i, m c ). Making 
an approximation m s — m d m s , Eq.(27) can be simplified as 

(fB s /fB d )/(foJfo d ) = 1 + - a 2 ^. (28) 
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The small current quark mass allows us to evaluate precisely fx and /s at m q m 0. Along 
the same line as in Ref . , it is straightforward to get the sum rules for f K and /s 



For simplicity, we do not give the final expressions for these two coefficients. To calculate these 



relation between the threshold u c and the light quark mass m q is unknown. A resonable and 
simple assumption for this relation, according to the definition of the u in HQET, is u c (m q ) = 
lv c (0) + Am q with a coefficient A varying between 0.0 and 2.0. In fact, this linear relation has 
successfully be applied to QCD sum rule calculations in Ref. p. Using the same energy scale 
jj, = lGeV and condensate parameters as in Ref. 0, we can carry out the sum rule calculations 
for fa and /s, and further obtain our desired results for a\ and a 2 . In the manipulation process, 
we find a stable window T = 0.7 ~ 1.0 GeV to exist for the threshold u c = 2.0 ± 0.3 GeV. Our 
numerical result shows that a\ and a 2 are insensitive to variation of A from 0.0 to 2.0, because 
of large cancellation effects between the first two terms in the expressions for a x and a 2 . For 
instance, when taking A = 0.0, 0.5, 1.0, 1.5 and 2.0, we have a x = -1.188, -1.201, -1.208, -1.213 
and -1.215, and a 2 = -0.905, -0.916, -0.921, -0.924 and -0.925, respectively (an error of 4% in 
each exists due to the uncertainty of the comtinuum threshold u c ). In addition, it is worthwhile 
to emphasize that the uncertainties of the quark masses produce only the corrections of about 
2% to ai and a 2 . The two coefficients ai and a 2 in case of A=l, as the functions of the Borel 




(29) 



and 




(30) 



two coefficients at m q ~ 0, the formulas^ for differential calculus can be used. However, an exact 
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parameter T, are plotted in Fig.l. 

To summarize, in this letter, we discuss the double ratio (fB s /fB d )/(fD s /fD d ) by applying 
QCD sum rules in HQET. We find that the correction term ai(m s /mb — m s /m c ) for the double 
ratio in Ref.|| and Ref.jlJ is modified to ax(m s /mb) — a2{m s /m c ) slightly, with a\ ~ — 1.2 and 
a 2 ~ —0.9. Our expression for the double ratio differs from those of Grinstein and Oakes in two 
aspects: (i) the double ratio in our approach is determined not only by quark mass terms but 
also by vacuum condensates, while in Ref.|| and Ref.0] the SU(3) breaking effect in the double 
ratio results only from the light quark masses; (ii) as a result in HQET, a\ and a 2 depend weakly 
on the heavy quark masses and m c , respectively, via a s (mb) and a s (m c ). This leads to the 
small difference between a\ and a 2 . Considering uncertainties of the heavy meson and the quark 
masses, we can conclude that the double ratio is a little greater than 1. This is also consistent 
with Oakes's prediction. Our conclusion will be beneficial to precision determination of fsj fB d 
from the data on f Ds J f Dd in future r — c factory, and further has us get a good understanding 
of the B — B mixing phenomena. 

We would like to thank Professor R.J. Oakes for helpful discussions and introducing Ref.[[| 
to us. This work is in part supported by the National Science Foundation of China. 



9 



References 

[1] C.Bernard, J.Labrenz, and A.Soni, in Proceedings of Lattice '93, Amsterdam 1992 
(P.vanBaal and J.Smit,eds.), Nucl.Phys.B (Proc.Suppl.) 30,465 (1993); Phys.Rev.D49, 
2536(1994). 

[2] C.A.Dominguez, in Proceedings of the Third Workshop on the Tau-Charm Factory, Mar- 
bella, Spain, June, 1993. 

[3] B.Grinstein, Phys.Rev.Lett.71, 3067(1993). 

[4] R.J.Oakes, Phys.Rev.Lett. 73,381(1994). 

[5] Tao Huang and Chuan-Wang Luo, Phys.Rev.D53,5042(1996). 

[6] H.Georgi, Phys.Lett. B240,447(1990); J.G.K6rner and G.Thompson, Phys. Lett. B 
264,185(1991); T.Mannel, W.Roberts, and Z.Ryzak, Nucl. Phys. B 368, 204(1992). 

[7] M.Neubert, Phys.Rev.D 45,2451(1992);46,1076(1992). 

[8] M.A.Shifman, A.I.Vainshtein and V.I.Zakharov, Nucl. Phys. B147, 385 (1979); B147, 
448(1979). 



10 



